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We present the results of a search for gravitational waves associated with 223 y - ray bursts (GRBs) 
detected by the Interplanetary Network (IPN) in 2005-2010 during LIGO’s fifth and sixth science runs 
and Virgo’s first, second, and third science runs. The IPN satellites provide accurate times of the bursts 
and sky localizations that vary significantly from degree scale to hundreds of square degrees. We search 
for both a well-modeled binary coalescence signal, the favored progenitor model for short GRBs, and 
for generic, unmodeled gravitational wave bursts. Both searches use the event time and sky localization 
to improve the gravitational wave search sensitivity as compared to corresponding all-time, all-sky 
searches. We find no evidence of a gravitational wave signal associated with any of the IPN GRBs in 
the sample, nor do we find evidence for a population of weak gravitational wave signals associated with 
the GRBs. For all IPN-detected GRBs, for which a sufficient duration of quality gravitational wave data 
are available, we place lower bounds on the distance to the source in accordance with an optimistic 
assumption of gravitational wave emission energy of 1 0 _2 M Q c 2 at 150 Hz, and find a median of 
13 Mpc. For the 27 short-hard GRBs we place 90% confidence exclusion distances to two source 
models: a binary neutron star coalescence, with a median distance of 12 Mpc, or the coalescence of a 
neutron star and black hole, with a median distance of 22 Mpc. Finally, we combine this search with 
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previously published results to provide a population statement for GRB searches in first-generation 
LIGO and Virgo gravitational wave detectors and a resulting examination of prospects for the advanced 
gravitational wave detectors. 

DOI: 10.11 03/Phy sRevLett. 1 1 3 .0 1 1 1 02 PACS numbers: 04.80.Nn, 07.05.Kf, 95.85.Sz, 98.70.Rz 


y-ray bursts (GRBs) are amongst the most energetic 
electromagnetic astrophysical phenomena. They fall into 
two commonly accepted groups depending on their dura- 
tion and spectral hardness [1,2] and are referred to as 
“short” or “long.” Short GRBs (duration less than 2 s, hard 
spectra) are believed to be produced by the mergers of 
either double neutron star or neutron-star-black-hole bina- 
ries [1], and the recent observation of a kilonova associated 
with GRB130603B [3,4] lends support to this hypothesis. 
Such compact binary coalescences generate strong gravi- 
tational waves (GWs) in the sensitive frequency band of 
Earth-based gravitational wave detectors [5,6]. The detec- 
tion of gravitational waves associated with a short GRB 
would provide direct evidence that the progenitor is indeed a 
coalescing compact binary as well as distinguish between a 
double neutron star (NSNS) and neutron-star-black-hole 
(NSBH) progenitor. Long GRBs (duration greater than 2 s, 
soft spectra) models are mostly related to the collapse of 
rapidly rotating massive stars. The extreme conditions 
encountered in such objects may make the system suscep- 
tible to a variety of rotational instabilities that may emit up to 
10 _2 M o c 2 through GW radiation [7,8]. 

Between 2005 and 2010, the first generation of large- 
scale interferometric gravitational wave detectors — LIGO, 
Virgo, and GEO — were operating at, or close to, their 
design sensitivities. During these runs, the detectors had a 
sensitivity out to tens of Mpc for binary mergers [9-12] and 
other transients emitting 10 _2 A/ Q c 2 in gravitational waves 
[13-15]. 

Although it is expected that most GRB progenitors 
observed by y - ray detectors will be at distances too large 
for the resulting gravitational wave signals to be detectable 
by initial LIGO and Virgo [16,17], it is possible that in 
the GRB data set under study that one might be within the 
range of the detectors. For example, the smallest observed 
redshift to date of an optical GRB afterglow is z — 0.0085 
(—36 Mpc) for GRB 980425 [18-20]; this would be within 
the LIGO- Virgo detectable range for some progenitor 
models. Although GRB 980425 is a long duration soft 
spectrum GRB, observations seem to suggest that, on 
average, short-duration GRBs tend to have smaller redshifts 
than long GRB s [2 1 ,22] . We therefore search for evidence of 
a gravitational wave signal associated with any observed 
short or long GRB for which there is a sufficient duration of 
high-quality data in at least two detectors. By making use of 
the known time and sky location of the observed GRB, it is 
possible to significantly reduce the parameters of the search 
and consequently improve the sensitivity over an all-sky, all- 
time search of the data. Several searches for gravitational 


waves associated with y - ray bursts have been performed 
in the past using data from both LIGO and Virgo detectors 
[23-27]. Indeed, the data from LIGO’s fifth and sixth 
science runs (S5 and S6) and Virgo’s first through third 
science runs (VSR1-3) were analyzed to search for gravi- 
tational wave signals associated with both short and long 
GRBs observed with the Swift BAT and Fermi GBM and 
L AT detectors [28-30] . No evidence for a gravitational wave 
signal was found in these searches. 

The Interplanetary Network (IPN) [31] is a group of 
satellites orbiting the Earth, Mars, and Mercury and 
operating, among other equipment on board, y - ray detec- 
tors. The IPN, in its current configuration, acts as a quasi- 
all-sky and full-time y - ray burst detector. Thus, the IPN 
provides an additional population of GRBs which may not 
be observed solely by Swift or Fermi and tends to detect 
brighter (therefore, on average, closer) bursts, which are 
relevant for gravitational wave searches. The IPN provides 
accurate GRB times as well as sky locations determined by 
triangulation between the satellites in the network. 
Depending upon the satellites involved, the localization 
can vary from less than one square degree to over one 
thousand. Two short duration GRBs, GRB 051103 [32,33] 
and GRB 070201 [34], were localized by the IPN with error 
boxes overlapping the M81 galaxy at 3.6 Mpc and the 
Andromeda galaxy (M31) at 770 kpc, respectively. The 
gravitational wave data around the times of these GRBs 
were analyzed in [35] and [36], respectively. The non- 
detection of associated gravitational waves ruled out the 
progenitor object being a binary merger in M81 or M31 
with over 99% confidence. 

In this Letter we present the results of a targeted search 
for gravitational waves around the burst trigger times of 
223 additional y - ray bursts, including 27 short GRBs, 
localized by the Interplanetary Network (IPN) during both 
LIGO’s fifth and Virgo’s first science runs (S5 and VSR1) 
and LIGO’s sixth and Virgo’s second and third science runs 
(S6 and VSR2-3). The search for gravitational wave bursts 
(GWBs) is performed on all the GRBs, short or long, for 
which we have good-quality gravitational wave data, 
regardless of the localization error box size. In addition, 
a search for a binary merger signal is performed for the 27 
short hard GRBs for which there was sufficiently good sky 
localization to make the search tractable. 

We find no evidence for a gravitational wave candidate 
associated with any of the IPN GRBs in this sample, and 
statistical analyses of the GRB sample rule out the presence 
of a collective signature of weak gravitational waves 
associated with the GRB population. We place lower 
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bounds on the distance to the progenitor for each GRB, and 
constrain the fraction of the observed GRB population at 
low redshifts. Additionally, we combine the results pre- 
sented here with those from previous analyses to provide a 
comprehensive limit from all of the GRB searches. Using 
this, we extrapolate to the future advanced detector network 
and show that the observation of a gravitational wave signal 
associated with a GRB is possible, but by no means 
guaranteed and the continued all sky coverage provided 
by the IPN will increase these chances. 

GRB sample . — The Interplanetary Network (IPN) is a 
group of spacecraft equipped with y - ray detectors used to 
localize y - ray bursts (GRBs) and soft y repeaters (SGRs, or 
magnetars). The IPN has contributed burst data to LIGO 
since the initial engineering runs in 2001 . At the time of this 
combined search, nine spacecraft contribute their data: 
Wind, Mars Odyssey, MESSENGER, INTEGRAL, RHESSI, 
Swift, Suzaku, AGILE and Fermi. 

The astronomical locations of GRBs are determined by 
comparing the relative arrival times of the signal at the 
spacecraft. The precision is inversely proportional to the 
spacecraft separations, among other things, so that 
the localization accuracy of a network with baselines of 
thousands of light-seconds can be equal to or superior to 
that of any other technique. A description of the error box 
construction process can be found in [37] and, specific to 
our search, in [38]. The light curves, energy spectra, and 
localizations of all the bursts in our sample were examined 
to eliminate the possibility of contamination by magnetar 
bursts or solar flares. None of these events have been 
followed up by x-ray or optical telescopes, so no informa- 
tion on afterglows or possible host galaxies and associated 
redshifts is available. The full list of GRBs and their 
parameters can be found in the Supplemental Material [39] 
associated with the Letter. 

The classification of GRBs into long or short typically 
uses the T 90 duration for a given detector. This is defined as 
the time interval over which 90% of the total background- 
subtracted photon counts are accumulated. However, this 
time depends on the energy range to which the detector is 
sensitive and may therefore vary across the satellites 
observing the bursts. Since the IPN bursts are observed 
by a set of different detectors with different sensitivities, to 
quote a single T 9() for them could be misleading. Even 
when a single detector measures this time, it is possible to 
get different numbers for the same burst depending on the 
arrival angle, which affects the sensitivity as a function of 
energy. In this analysis, where possible, we have used the 
classification provided by [40], based on observations with 
Konus-Wind. We note that the set of short bursts observed 
by Konus is split into type I (likely merger scenario), 
possibly with extended emission, and type II (collapsar). 
For the modeled search, we only analyzed the type I bursts. 
For bursts not observed by Konus, the T 90 observed by 
Suzaku was used and in cases where this was not available, 


a by-eye estimate of duration from another mission with 
good sensitivity (such as Swift or INTEGRAL ) was used. In 
these cases, any burst with a T 90 under two seconds was 
classified as short. The unmodeled search uses many 
minutes of data around the given time of the burst, whereas 
the modeled search uses a small six-second window around 
“short” bursts to search for the binary coalescence, due to the 
predicted time difference between merger and GRB emis- 
sion. It is therefore important that the time of arrival at the 
Earth be calculated as accurately as possible for these bursts. 
When the burst is observed by a near-Earth satellite, this is 
straightforward and approximated to the spacecraft time. 
However, for some poorly localized GRBs only observed by 
distant satellites, the uncertainty in the Earth-crossing time 
can be up to approximately five seconds. 

Only the GRBs that occurred when two or more of the 
LIGO and Virgo detectors were operating in a stable 
configuration are analyzed. Gravitational-wave data seg- 
ments that are flagged as being of poor quality are excluded 
from the analysis. Thus, although the IPN observed over 
600 bursts during the period of interest, only 223 could be 
analyzed, of which 27 are classified as short. This never- 
theless constitutes the largest GRB sample used to date in 
such a study. 

Gravitational wave detectors . — In this Letter, we discuss 
results obtained from analyzing data collected by the initial 
LIGO and Virgo gravitational wave detectors. There were 
three LIGO detectors: a 4 km and a 2 km detector, both at 
Hanford, Washington, (referred to as HI and H2, respec- 
tively) and another 4 km detector in Livingston, Louisiana 
(LI). The Virgo detector is a 3 km detector located in 
Cascina, Italy. Details of these detectors can be found in 
[41] and [42]. From 2005 to 2010, these detectors were 
operating at or near design sensitivity. The fifth LIGO 
science run (S5) took place from November 2005 to August 
2007 and the sixth science run (S6) ran from June 2009 to 
October 2010. The second Hanford detector (H2) was not 
operational during S6. Virgo operated three distinct science 
runs during this time: VSR1 ran from May to October 
2007, VSR2 ran from July 2009 to January 2010, and 
VSR3 ran from August to October 2010. 

The existence of multiple, widely separated detectors 
around the globe aids our ability to localize signals in blind, 
all-sky searches. Additionally, noise artifacts in the detectors 
caused by terrestrial disturbances are likely to be uncorre- 
lated between the detectors, aiding the ability to distinguish 
true signals from the noise background. The Hanford and 
Livingston detectors are separated by 3000 km, which 
corresponds to 10 ms travel time for light or gravitational 
waves. The Hanford and Virgo observatories are separated 
by 27 ms and the Livingston and Virgo observatories 
by 26 ms. 

Search methodology . — The methods used to search for 
binary merger signals and GWBs are largely the same as for 
the previous analysis described in [30]. The one major 
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change is the necessity to search the variably shaped, 
irregularly sized error regions provided by the IPN local- 
izations. We begin by describing how that is done and then 
provide a brief review of the remainder of the analysis 
details, referring the reader to [30] for more details. 

Covering the error boxes: The analysis of the gravita- 
tional wave data depends on the assumed sky direction to 
the GRB, since the data must be time-shifted according to 
the expected time-of-arrival at each detector. This is also 
weighted by the response of each GW detector to the 
assumed sky direction. Each GRB localization error box 
corresponds to a 3o region determined by the intersection 
of the IPN timing annuli, with a construction process 
described in detail in [38]. However, most of these IPN 
error regions are larger than the directional resolution of the 
GW detector network. Thus, to maximize the likelihood of 
finding a gravitational wave signal associated with the 
GRB, we perform a discrete search across the entire IPN 
error box by populating it with a grid of points and 
repeating the search at each of these points. The GW 
detectors have a timing resolution of ~0.5 ms, correspond- 
ing to a spatial resolution of a few degrees. Therefore, we 
chose search point separations of approximately 3.6 
degrees when only LIGO interferometers’ data were used 
in the search and 1.8 degrees when the more widely 
separated Virgo detector was also used. The probability 
distributions of search points over the error boxes were 
chosen Gaussian for long bursts (to assure that there are 
proportionally more points for those positions with larger 
probabilities of containing a signal) and uniform for short 
bursts (no assumptions made in terms of signal origin 
within the error box). Short GRBs which could not be 
localized to better than a few hundred square degrees were 
not analyzed in the modeled search due to high computa- 
tional requirements and the negligible increase in sensi- 
tivity rendered by a targeted search. We use simulated 
signals to determine the sensitivity of our searches. We 
distribute these over the IPN error boxes, with the density 
of simulations weighted according to the estimated source 
position probability distribution. This assures that there are 
proportionally more simulations for those positions with 
larger probabilities of being the true GRB signal location. 

During the S5 run, two LIGO detectors, HI and H2, were 
operational at the Hanford site. There are a number of 
bursts for which the only available gravitational wave data 
are from these two detectors. Since the detectors are co- 
located and co-aligned, it is not possible to distinguish 
different sky locations based on the observed gravitational 
wave signal. In this case, there is no need to perform an 
explicit search over the error box and, for the binary search, 
all such GRBs can be analyzed. 

Search for GWs from a compact binary progenitor: For 
the binary merger progenitor model, it is believed that the 
delay between the merger and the emission of y rays will be 
small [see discussion in [30]]. We therefore search for 


binary coalescences with a merger time between 5 s prior to 
the GRB and 1 s afterwards. This is wide enough to allow 
for potential precursors, some uncertainties in the emission 
model and in the arrival time of the electromagnetic signal 
at the IPN spacecraft, as well as for the differences in 
sensitivity of the IPN detectors. In addition, we require a 
minimum of 40 minutes of data available around the time of 
the GRB. This ensures both that the detectors were 
operating stably at the time as well as providing a set of 
comparable data which can be used to estimate the back- 
ground of noise events. 

For the short GRBs, the data streams from the opera- 
tional detectors are combined coherently and searched by 
matched-filtering against a bank of binary merger gravita- 
tional waveforms, as described in [43]. The gravitational 
waveform emitted by a binary system depends on the 
masses and spins of the NS and its companion (either NS or 
BH), as well as on the distance to the source, its sky 
position, inclination angle, and the polarization angle of the 
orbital axis. As described above, we tile the sky region with 
a fixed set of points to search, and we similarly tile the mass 
space [44] to provide sensitivity to binaries with component 
masses greater than 1 M Q with an upper limit of 3 M e for 
any neutron star and a total mass of 25 M 0 . The remaining 
parameters are handled by maximizing the likelihood 
analytically over these dimensions. In addition to matched 
filtering, the analysis utilizes a number of signal consis- 
tency tests to reject nonstationary, transient noise “glitches” 
in the GW detectors’ data [43]. 

Search for gravitational wave bursts: The procedure 
used to search for generic short-duration (< 1 s) GW bursts 
follows that used in previous GRB analyses [30,45]. All 
GRBs are treated identically regardless of their classifica- 
tion. We search for a GW event between 600 s prior to the 
GRB trigger time and 60 s after or the T 90 time, whichever 
is greater. This timescale allows us to take into account 
almost all of the possible scenarios for a GW emission 
associated with the GRB; see [30] for more details. Since 
the GWB search requires more data around the GRB than 
the binary merger search, there are GRBs for which we can 
perform the merger search but not the search for unmodeled 
transients. The data within a ±1.5 h window around the 
GRB is used to estimate the background of noise events in 
the data. The search for a GWB between 60 and 500 Hz is 
performed by the x-pipeline algorithm [46,47]. Candidate 
events are then ranked based on their energy and a number 
of signal consistency tests are applied to reduce the effect of 
nonstationary noise seen in the GW detectors. Events 
occurring at times of known instrumental problems or 
environmental disturbances are discarded. 

Significance of results: The significance of any candidate 
gravitational wave event is estimated by using the data 
surrounding the GRB time. Specifically, we divide the 
surrounding data into a large number of blocks of identical 
length to the search region and calculate the false alarm 
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probability (FAP), or p value of the event by counting the 
fraction of off-source trials with an event louder than the 
one observed. Where necessary, we artificially time-shift 
the data from the detectors by several seconds to generate 
more off-source trials that can be used to estimate the 
background. 

In addition to a single, significant event, it’s possible that 
the data could contain several weak signals. In order to test 
whether this is the case, we use a weighted distribution of 
observed p values and see whether it is consistent with the 
expected, uniform distribution of noise. We also use a 
weighted binomial test to more quantitatively assess con- 
sistency with the no-signal hypothesis. This test looks for 
deviations from the null hypothesis in the 5% tail of lowest 
p values weighted by the prior probability of detection 
(estimated from the GW search sensitivity). This combi- 
nation allows us to give more weight to those GRBs for 
which the gravitational wave network was most sensitive, 
which are the ones we are most likely to detect. The result 
of the weighted binomial test is compared with the 
distribution obtained from simulated results with p values 
uniformly distributed in [0, 1] to evaluate the population 
significance. The test is described in greater detail in the 
appendix of [30] and in [47]. 

Results. — A search for gravitational waves has been 
performed for a total of 223 GRBs. Of these, 27 were 
classified as short GRBs with a likely binary coalescence 
progenitor. The gravitational wave data at the time of these 
GRBs has been searched for evidence of the coalescence 
waveforms. For 221 bursts, including 25 of the 27 short 
bursts we have performed a search for generic gravitational 
wave transients in the IPN error box around the time of the 
GRB. A full list of all GRBs analysed is available in the 
supplementary material, where we provide two tables 
listing the short and long GRBs analyzed in this search. 

Modeled coalescing binary search results: For each 
GRB, we estimate the p values of the gravitational wave 
candidate, by comparing with the background trials. The 
distribution of observed p values is shown in Fig. 1. No 
significant candidates found. The p value is estimated from 
the background. However, for a number of GRBs, particu- 
larly those observed in the colocated Hanford detectors, the 
search yields no candidate gravitational wave events after 
background rejection cuts. For such GRBs, when no event 
is observed, we cannot quote an exact p value but only a 
range bounded below by the fraction of background trials 
with an event and above by 1. The result of the weighted 
binomial population detection test yields a background 
probability of « 98%, strongly favoring the no-signal 
hypothesis. In conclusion, no noteworthy individual events 
were found by this search, nor evidence for a collective 
population of weak gravitational wave signals. 

Unmodeled GW burst results: The unmodeled gravita- 
tional wave burst pipeline analyzed both short and long 
GRBs. The distribution of p values for each of the 221 IPN 



p value 


FIG. 1 (color online). Cumulative distribution of p values from 
the analysis of 27 short-duration IPN GRBs for evidence of a 
binary merger gravitational wave signal. The expected distribu- 
tion under the no-signal hypothesis is indicated by the dashed 
line. For GRBs with no event in the on-source region, we provide 
upper bounds on the p value equal to 1. 

GRBs analyzed is shown in Fig. 2. The binomial test yields 
a background probability of 68%, consistent with the null 
hypothesis. The smallest p value, 0.5%, came from 
GRB060203B, which is statistically consistent with a 
no-signal hypothesis given the number analyzed. 

Astrophysical interpretation . — Given that no significant 
event was found in our analyses, we place limits on GW 



FIG. 2 (color online). Cumulative distribution of p values from 
the analysis of 221 IPN GRBs for evidence of a gravitational 
wave transient associated to the burst. The expected distribution 
under the no-signal hypothesis is indicated by the dashed line. 
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FIG. 3 (color online). Histograms across the sample of short 
IPN GRBs of the distance exclusions at the 90% confidence level 
for NS-NS and NS-BH systems. 

emission based on both binary merger (for short GRBs) and 
generic gravitational wave burst (for all GRBs) signal 
models, and assess the potential of a similar search with 
second-generation gravitational wave detectors around 
2015-2020. 

Distance exclusion: For a given signal morphology, the 
gravitational wave analysis is efficient in recovering signals 
up to a certain distance limit that depends on the sensitivity 
of the detectors at the time of the search. We quote a 90% 
confidence level lower limit on the distance to each GRB 
progenitor: that is, the distance at which we recover 90% of 
simulated signals. The quoted exclusion distances are 
marginalized over systematic errors that are inherent in 
this analysis: errors introduced by the mismatch of a true 
GW signal and the waveforms used in the simulations [9] 
and amplitude errors from the calibration of the detector 
data [48]. 

For the short GRBs, we calculate a distance exclusion for 
both two neutron stars (NSNS) and a neutron star with a 
black hole (NSBH). In both cases, we assume a jet half- 
opening angle of 30°, and assume that the GRB is emitted 
in the direction of the binary’s total angular momentum. 
The median exclusion distance for NSNS is 12 Mpc and for 
NSBH is 22 Mpc. A histogram of their values is shown in 
Fig. 3. The NS masses are chosen from a Gaussian 
distribution centered at 1.4M Q [49,50] with a width of 
0.2 Mq for the NSNS case, and a broader spread of 0.4 M Q 
for the NSBH systems, to account for larger uncertainties 
given the lack of observations for such systems. The BH 
masses are Gaussian distributed with a mean of 1 04/ Q and 
a width of 6 M Q . The BH mass is restricted such that the 
total mass of the system is less than 25 M Q . For masses 
greater than this distribution, the NS would be swallowed 
without disruption by the BH, no massive torus would 
form, and no GRB would be produced [51-53]. 
The dimensionless NS spins are drawn uniformly over 


io s iff' io° to’ 10 * 

exclusion distance (Mpc) 

FIG. 4 (color online). Histograms across the sample of IPN 
GRBs of the distance exclusions at the 90% confidence level for 
circularly polarized sine-Gaussian GW burst models at 150 
and 300 Hz. We assume an optimistic standard siren GW 
emission of I? GW = 10 _2 M o c 2 . 

[0, 0.4], and the BH spins are drawn uniformly over 
[0, 0.98) with tilt angle < 60°. 

For the GWB search, no specific waveform model is 
assumed. Consequently, we use generic signal morpholo- 
gies to give an idea of the search sensitivity. Specifically, 
we use circularly polarized sine Gaussians with central 
emission frequencies of 150 and 300 Hz. We assume a jet 
opening angle of 5°, which is appropriate for long GRBs. 
We also assume a total GW emission of 10 _2 M o c 2 ; this 
corresponds to the most optimistic models for gravitational 
wave emission from long GRBs [30]. The median exclu- 
sion distance is 13.0 Mpc at 150 Hz and 4.9 Mpc at 300 Hz, 
and histograms of the distributions are given in Fig. 4. 



FIG. 5 (color online). Cumulative distribution of p values from 
the analysis of 508 GRBs. The expected distribution under the 
no-signal hypothesis is indicated by the dashed line. 
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FIG. 6 (color online). Cumulative distribution of p values for 
all 69 analyzed short GRBs. The expected distribution under the 
no-signal hypothesis is indicated by the dashed line. For GRBs 
with no event in the on-source region, we provide upper bounds 
on the p value of 1. 

Population exclusion from all GRBs analyzed: Here we 
present the combination of all S5-6 and VSR1-3 searches 
for coincident GRB and GW signals from this Letter, the S5 
GRB search [28,29] and the recent S6 and VSR2-3 search 
[30]. Algorithms which were used for the first S5 paper 
were adjusted and reviewed to make sure all analyses were 
comparable. In total, 508 GRBs were analyzed with the 
burst search and 69 short GRBs were analyzed for a 
compact binary coalescence GW signal. None of the 
separate searches showed evidence of a population of weak 
events. Figures 5 and 6 show the distribution of p values of 
the full set of GRBs. The weighted binomial test applied to 
the full population of GRBs confirms that the observed 
distributions are consistent with the null hypothesis (no 
signal being observed). 

Next, we use the full population of GRBs to place 
exclusions on GRB populations. To do this, we use a simple 
population model, where all GRB progenitors have the 
same GW emission (standard sirens), and perform exclu- 
sion on cumulative distance distributions. We parameterize 
the distance distribution with two components: a fraction 
F of GRBs distributed with a constant comoving density 
rate up to a luminosity distance R, and a fraction 1 — F at 
effectively infinite distance. This simple model yields a 
parameterization of astrophysical GRB distance distribu- 
tion models that predict a uniform local rate density and a 
more complex dependence at redshift > 0.1, as the large- 
redshift part of the distribution is well beyond the sensi- 
tivity of current GW detectors. The exclusion is then 
performed in the ( F , R) plane. For details of this method, 
see Appendix B of [30]. 





FIG. 7 (color online). Cumulative redshift distribution F(R) 
exclusion from the analysis of 508 GRBs with the GW burst 
search. We exclude at 90% confidence level cumulative distance 
distributions which pass through the region above the black solid 
curve. We assume a standard siren sine-Gaussian GW burst at 
150 Hz with an energy of £ G w = 10 _2 M o c 2 . We extrapolate this 
exclusion to Advanced LIGO and Virgo detectors assuming a 
factor 10 improvement in sensitivity and a factor 2 increase in 
number of GRB triggers analyzed. The blue dashed curve is the 
extrapolation assuming the same standard siren energy of Ii GW = 
10 _2 MqC 2 and the green (gray) dashed curve assuming a less 
optimistic standard siren energy of Ii GW = 10 _4 M q c 2 [54,55]. 
For reference, the red staircase curve shows the cumulative 
distribution of measured redshifts for Swift GRBs [56,57]. 

In Fig. 7 we show the exclusion for GW bursts, using as a 
reference signal a 150 Hz sine-Gaussian signal, with an 
energy in gravitational waves of 10 -2 A/ 0 c 2 . In addition, 
we plot the redshift distribution of GRBs as observed by 
Swift. The exclusion at low redshift is dictated by the 
number of GRBs analyzed and at high redshift by the 
typical sensitive range of the search. These exclusions 
assume 100% purity of the GRB sample. In Figure 8, we 
show the exclusion for the NSNS and NSBH sources. In 
neither case does the exclusion line come close to the 
observed population redshift, indicating that it would have 
been unlikely to observe an event in this analysis. Indeed, 
an analysis of all IPN bursts shows that their average 
redshift is 1.7, and that it detects short bursts with good 
efficiency up to a redshift of about 0.45. 

The advanced LIGO and Virgo detectors are nearing 
completion. They are expected to start taking data in 2015 
and reach their design sensitivity towards the end of the 
decade [58]. We can use the results obtained here to 
extrapolate and predict what might be expected with the 
advanced detectors. The S5-6 and VSR1-3 runs comprised 
a total of around 21 months of two (or more) detector duty 
cycle. Over that period, the detectors’ reach varied by 
approximately a factor of 4, from a 5 Mpc sensitive distance 
to NSNS sources for H2 in early S5 to 20 Mpc for HI and 
LI by the end of S6. Similarly, the current scenario [58] 
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FIG. 8 (color online). Cumulative distance exclusion for 69 
analyzed short GRBs for both a NSNS and a NSBH progenitor 
model. The exclusion distance is given for this test then 
extrapolated by a factor of two in number and ten in sensitivity 
for the advanced detector era expectations. For reference, the red 
staircase curve shows the cumulative distribution of measured 
redshifts for Swift GRBs [56,57]. 

calls for around 18 months of science runs of ever 
increasing sensitivity during the commissioning phase, 
prior to extended running at design sensitivity which will 
provide a reach roughly ten times what was achieved in S6 
and VSR2-3. 

To approximate the expected advanced detector results, 
we scale the exclusion distances obtained here by a factor 
of ten and also increase by a factor of two the number of 
observed GRBs to account for the increased run time of a 
few years. These extrapolated curves are also shown on 
Figs. 7 and 8. We see that for both generic burst signals and 
binary mergers, the exclusion curves are now comparable 
with the observed redshifts. In the optimistic scenario 
where every GRB emits 10 _2 M o c 2 of energy in gravita- 
tional waves, we will expect to see several such signals, or 
alternatively be able to exclude this scenario. However, 
since the 10~ 4 M e c 2 predicted exclusion is to the left of the 
observed redshifts, we will be unlikely to observe a 
gravitational wave associated to a GRB if this is the typical 
amplitude. For binary mergers, the NSNS and NSBH 
extrapolations bracket the line of observed redshifts, 
indicating that we might expect to see one (or fewer) 
NSNS associated with a GRB, but might expect several if 
NSBH are the progenitors. (These extrapolations are 
broadly comparable to those obtained using only the S6 
and VSR2-3 [30]. They are, however, slightly more 
pessimistic as they use a more realistic estimate of the 
evolution of detector sensitivity, as described above.) 

Of course, all of the above relies heavily on the 
continued operation of all sky sensitivity to GRBs, and 
reasonable localization ability. Swift and Fermi will con- 
tinue to run, and SVOM is expected in the advanced 


detector timeline, but it’s clear that the addition of the 
all-sky, full-time IPN will aid these searches as well. 

Conclusion . — Wc have performed a search for gravita- 
tional waves coincident with 223 y-ray bursts localized by 
the Interplanetary Network between 2005 and 2010. These 
GRBs were detected by the IPN during LIGO’s fifth and 
sixth science runs and Virgo’s first, second, and third 
science runs. Of these, we analyzed the 27 short GRBs with 
a focused search that looked for gravitational wave signals 
from the merger of neutron-star-neutron-star or neutron- 
star-black-hole binaries, as most likely expected for short 
GRBs. We also performed an unmodeled burst search over 
221 GRBs, both short and long. No gravitational wave was 
detected in coincidence with a GRB, and lower limits on 
the distance were set for each GRB for various gravitational 
wave emission models. 

Finally, we have combined these results with those of 
previous analyses for GRBs during S5-6 and VSR1-3 to 
provide a comprehensive statement on gravitational wave 
emission by GRBs from the first-generation LIGO and 
Virgo detectors. We also extrapolate this exclusion distance 
for the advanced detector era, where we assume a factor of 
2 increase in GRBs and a factor of 10 improvement in 
sensitivity. This shows that the advanced detector era will 
begin to exclude certain expected models and possibly 
make coincident gravitational wave detections with GRBs. 
These results and prospects for the advanced detector era 
demonstrate the benefit of searches triggered by y-ray burst 
observatories with high sky coverage such as the 
Interplanetary Network. The continued operation of these 
satellites will be crucial in future coincident GRB and 
gravitational wave searches and increase the likelihood of a 
detection in the advanced detector era. 
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